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We use a modified directed path algorithm to study the finite temperature chiral singularities of 
two color lattice QCD with staggered fermions at strong coupling. Our lattice calculations are 
done at a fixed finite temperature in the broken phase with a variety of different quark masses. 
We find that to the lowest order, the behavior of our observables, namely the condensates and the 
decay constants are all consistent with the predictions of the low energy effective theory of our 
model. We further notice that in order to see consistency, the quark masses used in our lattice 
calculations need to be quite small than typically used. 
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1. INTRODUCTION 

To study the chiral singularities in lattice QCD is a challenge due to a variety of technical diffi- 
culties, for example, the lack of algorithms which work efficiently for the simulations with realistic 
pion mass and the logarithmic nature of the chiral singularities. Because of these difficulties, today 
most lattice simulations are performed using large quark masses and the results are then extrapo- 
lated to the physical regime by using chiral perturbation theory (ChPT). However, some attempts 
to connect lattice QCD with the predictions from ChPT do not give clear answers [Q, p[ ^| f|, 
Given the difficulties described above in understanding the chiral singularities of realistic QCD 
simulations, we focus on a simpler toy model, namely the two color lattice QCD at strong cou- 
plings. Further, since it is easier to detect power-like singularities as compared to logarithmic ones, 
we have done our lattice calculations at a fixed finite temperature in the broken phase. In order to 
study the finite temperature singularities of the model, we have developed an efficient algorithm 
based on the directed path algorithm With this algorithm, the singularities of our model can be 
studied quantitatively and the results will be present in the following.. 

2. The Model and the Symmetries 

The lattice action of two color lattice QCD (2CLQCD) model at zero baryon density which 
we have studied is given by 



Here m is the quark mass and the quark fields x(x) and %(x) are two components vector fields 
with Grassmann variables. The gauge fields U a (x) are elements of Si/ (2) group and live on the 
links between x and x + a where a = t, 1,2,3. The factor r a = 1 for a = 1,2,3 and r, = \/a t . 
By choosing an L t x L d lattice (periodic in all directions) we can study thermodynamics in the 
L — > oo limit at a fixed L t by defining T = \/{a t ) 2 as the parameter that represents the temperature 
of the model. The absence of the gauge action shows that we are in the strong gauge coupling 
limit. A detailed discussion of the symmetries of 2CLQCD can be found in [0, g]. Here we only 
review them briefly. When m = our model has a U (2) global symmetry. In particular, this U (2) 
global symmetry contains a Ug(l) baryon number symmetry and a.Uc(l) chrial symmetry. At low 
temperature, the U (2) global symmetry is spontaneously broken down to f/fl(l) which leads to the 
existence of three massless Goldstone particles. In the presence of non-zero quark masses, only the 
baryon number symmetry survives and remains unbroken 1 . Further, once massive quarks 

are present, the three Goldstone particles acquire masses. One we call "pion" and the other two, 
which are degenerate in mass are called "diquark baryons". 

3. Dimer Representation of the Model 

One of the computational advantages of the strong coupling limit is that in this limit it is possi- 
ble to rewrite the partition function as a sum over configurations containing gauge invariant objects, 

'This Ub{1) baryon number symmetry will survive until the chemical potential fl is switched on and reachs a critical 
value He. 



r a ria (x) x(x)U a (x)x(x + a) - %(x + a)Ul {x)%{x) + m%X (*) 



(2.1) 
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namely monomers, dimers and baryons [Q, [K]]. A simple calculation shows that the partition func- 
tion Z of the model can be rewritten as: 

/ - y rry — (I)[M*)+IM*)UrW*)+M*)U. (3.1) 
{a:} x i 2 

where the variables n x = 0, 1,2, ki(x) = 0, 1,2, bi{x) = 1,0, — 1, / = t, 1,2,3 associated with each 
lattice site x are defined as the number of monomers, dimer bonds and oriented baryon bonds 
connecting lattice sites x and x + 1 respectively. Further, the baryon bonds always form isolated 
self-avoiding closed loops and each lattice site x in the remaining configuration without baryon 



loops satisfies the constraint: + n x = 2. Note that the partition function in (3.1) is positive 

definite, hence a Monte-Carlo algorithm can in principle be designed to study this system. The 
algorithms used in our lattice simulations is an extension of the directed path algorithm discovered 
in We refer readers to [[jl]] and do not provide the details of this algorithm here. 

4. Observables 

A variety of observables can be measured with our algorithm. We will focus on the fol- 
lowing: (1) The chiral two point function and the chiral susceptibility Xc which are given by 

Gc{z,z') = (%(z)x{z)x(z')x{z , )\ and Xc = L z ' G c(z,z'). (2) The diquark two point function 

which is given by Gb{z,z') = l%i{z)%2(z)%z(£)%i(z!)\. (3) The helicity modulus associated with 

the 1/(1) chiral symmetry Y c = 5^1^=1,2,3 ((Lx^(x)) 2 ^j where A M (x) = s(x) [|& M (*)| +fy(x)] 
and V s = L 3 . (4) The helicity modulus associated with the f7(l) baryon number symmetry Yb = 

5^1^=1,2,3 \(I.xM x )) 2 J where = [b^x)]. 

In addition to these observables, the correlation lengths M % and Mb associated with the chiral 
and diquark two point functions are also measured by calculating the temporal and spatial decays 
of the related two point functions respectively. 

5. Results 

5.1 The Low Energy Effective Theory of the Model 

In order to study the finite temperature chiral singularities of the model, we have done exten- 
sive lattice simulations at a fixed finite temperature in the broken phase with a variety of different 
quark masses m. Further, since from our lattice calculations in the chiral limit, we found that the 
transition temperature Tq ~ 2.9285 [[TT|], hence we focus on the chiral singularities at T = 2.9, a 
temperature which is below and close to Tq- 

Before describing the details of our lattice results, we review the low energy effective theory of 
the model. For small quark masses, it is possible to argue that the effective action Sh of the model 
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Figure 1: (Left) This figure shows Yq and Yg as a function of box size L. The solid lines are the fits of Yq 
and Y B to formula (5.2). These fits shows that F% = 0.1434(8), Fg = 0.1251(7). (Right) This figure shows 
Xc as a function of box size L. The E can be computed by fitting %c to formula (5.3) (solid line) and is given 
by £ = 0.4146(8). 



is given by: 

Si, 



r i F 2 F 2 1 

/ d\l -j-id^S) ■ (d^S) + ^-{d^u) ■ {d^u) --LS0u-h\. (5.1) 

where S is an unit 3-vector, u is an unit 2-vector and h is given by h = 2ra(l,0,0) (g) (1,0). The F 2 , 
Fg and £ are the low energy constants of the effective theory 2 . Further, these low energy constants 
can be computated from our lattice calculations in the massless theory. For example, F 2 and Fg 



can be calculated by fitting Yq and Y% to their finite size scaling (FSS) formulae respectively Qll|]: 

Y c = F 2 + a'L- l +a"L- 2 , 



2 -F 2 + b'L- l +b"L- 2 . (5 - 2) 



where a', a", b' and b" are some constants. Further, using the extrapolated F 2 and Fg, £ can be 
obtained from the FSS formula of Xc'- 



T 2 
~6 



+ a\L. (5.3) 



where a\ is some constant depending on the low energy constants of the theory. By this method, 
at T = 2.9, we found that F 2 = 0.1434(8), Fj = 0.1251(7) and £ = 0.4146(8) (Figjl]). All the 
% 2 /DOF of fits in above and the following discussions are less than 1. On the other hand, these 
low energy constants can also be computed from our lattice simulations with massive quarks by 
performing the chiral extrapolation in the appropriate regime. In the follwing, we will discuss this 
procedure of chiral extrapolations in detail. 

5.2 The Pion Decay Constant 

In the regime where the condition LM K 3> 1 is satisfied (M n is the pion mass), one would 
expect that the observables such as the pion decay cosntant and condensate satisfy the expansion 



(O) = a Q + b Q x + c x 2 + d x 3 + 0(x 4 ). (5.4) 
2 In the following, we will call F K the pion decay constant and E the condensate 
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= 0.1621(18), a 2 = 8.5(5), b 2 = -25(3) f/ = 0.1445(8), a 2 = 17.5(11), b 2 = -104(13) 




x x 

Figure 2: (Left) This figure shows F% as a function of x. The solid line is the fit of F% with 0.06672 < x < 
0.1887 to the formula described in the text. This fit shows that F% — 0.1621 (18) which is not consistent with 
our e regime lattice calculations. (Right) The fit of F% with 0.0211 < x < 0.0731 shows F% = 0.1445(8) 
which matchs nicely with the prediction. The most left data points (red squares) in both figures are the 
extrapolated results described in the text and are added to both figures by hands. 



where x is a dimensionless parameter and is given by x = -^pp, oq is the same observable obtained 
from the calculations with zero mass (see above) and yjm is the power-like singularity arising 



due to the infrared pion physics [J12J. Further, from the effective lagrangian given by ( |5.1| ), one 



can argue that the coefficient b in ( p.4| ) vanishes for the pion decay constant [ |13| , |14j. Indeed, 
in the presence of massive quarks with mass m, if the pion decay constant F 2 is extrapolated by 
F% = lim^oo Yc, then we found that the x-dependent prediction F 2 + c\x 2 + dix 3 of the pion decay 
constant is captured nicely by these extrapolated F 2 . However, if we focus only on these F 2 with 
0.06672 < x < 0.1887, then the fit shows that F 2 = 0.1621(18) (left figure in Fig.|) which is 
obviously not consistent with F% = 0.1434(8) from our lattice calculations with massless quarks. 
On the other hand, if we restrict ourselves to the F 2 with 0.0211 < x < 0.0731, then the same 
analysis shows that F% = 0.1445(8) (right figure in Fig.|2|) which agrees nicely with its predicted 
value. If we fit F% to F% + b x x + qx 2 + dix 3 with 0.021 1 < x < 0.0731. Then the fit shows 
that F% = 0.1413(48) and b\ = 0.24(35). Again both of them are consistent with the predictions 
F% = 0.1434(8) and^j =0 

5.3 The Condensate 

In terms of the dimensionless parameter x, the condensate £ m obtained from our lattice simula- 
tions with massive quarks, is given by Z 2 = a2 +b2X + 0(x 2 ). Further, it can be shown that «2 = £ 2 
and b 2 = -(^)-i(^G^(O) + 4|G B (0)), where G,(0) = H- Usin § the F l f b and 

£ calculated in the massless theory, ci2 and b 2 are given by a 2 = 0.1719(10) and b 2 = 1.19(9) 
respectively. Again this prediction is captured nicely with the extrapolated Z 2 which is obtained 
through Lj n = lim^oo Xc/L 3 if we include data at small x. However, as what has been already 
shown in the case of pion decay constant, the fitting results obtained from £ 2 with larger x dose 
not show a consistent result (left figure in FigJ^J). Only with the £ 2 from the smaller x regime, the 
fit would show a satisfactory result, namely a 2 and b 2 calculated from the fit match nicely with our 
predictions (right figure in Fig .[3]). 
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— %PT 

L 2 = 0.1568(10) 
o, = 1.64(2) 








[iii 




2 2 

— E m = L + c, x 






— XPT 






L 2 = 0.1692(25) 






C, = 1.32(7) 
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Figure 3: (Left) This figure shows E 2 as a function of x. The blue solid line is the fit described in the text with 
0.05168 < x < 0.1 156 and the red solid line is the prediction of jPT The fit indicates E 2 = a 2 = 0.1568(10) 
and b 2 = 1 .64(2). Neither of them match our predictions. (Right) The fit of E 2 with 0.02497 < x < 0.05 168 
shows E 2 = 0.1692(25) and b 2 = 1.32(7) both of which match nicely with the predictions. Again, the blue 
solid line is the fit described in the text and the red solid line is the prediction of %PT. In both figures, the 
data points at m = are obtained from the e regime calculations (red squares) and the extrapolation results 
(blue squares). They are added to the figures by hands. 



6. Discussion and Conclusions 



We have used an efficient algorithm to study the chiral singularities of a two color lattice QCD 
model at a fixed finite temperature and have succeeded in understanding the low energy physics of 
the model in terms of a effective theory with a few low energy constants. In particular we found 
that the relation lim^oo lim m ^o O = lim m ^o hm^oo O is satisfied, where O is either the pion decay 
constant or the condensate. However in order to see consistency, small quark masses are needed. In 
addition to the pion decay constant and the condensate, we also have studied the pion and diquark 
masses which will be discussed elsewhere. 

This work is done in collaboration with Shailesh Chandrasekharan and we thank Prof. T. 
Mehen, Prof. R. Springer, Prof. U.-J. Wiese and Dr. B. Tiburzi for helpful comments. This work 
was supported in part by the Department of Energy (DOE) grant DE-FG-02-05ER41368. 
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